Macrospin model of incubation delay due to the field-like spin transfer torque 
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We show that the absence of pre-switching oscillations ("incubation delay") in magnetic tunnel 
junctions can be explained within the macrospin model by a sizable field-like component of the 
spin-transfer torque. It is further suggested that measurements of the voltage dependence of tunnel 
junction switching time in the presence of external easy axis magnetic fields can be used to determine 
the magnitude and voltage dependence of the field-like torque. 



A spin polarized electric current can transfer spin an- 
gular momentum to a magnetic material, generating a 
torque that can induce magnetization dynamics and even 
magnetization reversal [l|, |2j . While extensive measure- 
ments have tested the validity and limitations of the 
macrospin model with Slonczewski's spin transfer torque 
in metallic spin valves, recent experiments with magnetic 
tunnel junctions (MTJ's) 0,0, Hill 

have observed an ad- 
ditional "field-like" or "perpendicular" spin torque. The 
existence of a field-like torque was predicted for metal- 
lic spin valves @, H, @, but shown to be smaller than 
Slonkzewski's "parallel" torque [13, El- For magnetic 
tunnel junctions, however, it was predicted that both 
torques could have similar magnitudes and that the field- 
like torque would have a quadratic dependence on volt- 
age [H, EH- The observed field-like torques generally 
agree with theoretical predictions, but some controversies 
remain. For example,measurements in the frequency do- 
main at low voltages @,@| and measurements of switching 
currents at large voltages @ report contradictory signs 
of the field-like term. This suggests that further theoreti- 
cal analysis and experimental investigation are necessary 
to fully understand the origin and the functional form of 
this torque. 

Here we report that the field-like torque can explain 
the absence of the pre-switching oscillations ("incuba- 
tion delay") found by Devolder et al. [3 in MTJ's. This 
observation could not be described within a macrospin 
model of magnetization reversal based on Slonczewski's 
spin transfer torque alone. Such a model predicts that 
pumping of the ferromagnetic resonance mode produces 
increasing oscillations in the resistance before switch- 
ing [l5|, [l|| ■ However, by including the effects of a field- 
like spin torque term within a macrospin model we are 
able to reproduce the main features of the observed mag- 
netization reversal: (i) a slow regular change of the resis- 
tance without oscillations preceding the switching, (ii) 
decaying oscillations of the resistance after switching, 
and (iii) similarity between magnetization reversal curves 
shifted so as to align their switching times (in our case the 
switching times ts are distributed between and 10 ns). 
In addition, we propose time-domain experiments which 




FIG. 1: (a) Schematic of rectangular MTJ showing the free 
layer magnetization M at an arbitrary orientation, and the 
reference layer magnetization mp pinned along the easy axis. 
The polarity of the voltage V applied across the sample is 
shown, (b) Unit sphere describing the possible orientations 
of M. Unit vectors normal to the sample (e n ), and along the 
easy (en) and hard axis (ex) are shown. The orientation of M 
is described by 9, the angle between the free and pinned mag- 
netizations, and ip, the angle between e n and the projection 
of M onto the e„-ei plane. Points P and AP are the equilib- 
rium positions of M along the easy axis, (c), (d) Schematics 
of MTJ energy density landscape (c) before and (d) during a 
voltage pulse, for M in plane (ip — tt/2), where K is the easy 
axis anisotropy energy density, ksT is the thermal energy at 
temperature T, and v is the free layer volume. In (d) the 
field-like spin torque lowers the energy barrier allowing mag- 
netization switching via a combination of Slonczewski's spin 
torque and random fluctuations. 



could be used to measure the magnitude of the field-like 
spin torque term and its voltage dependence, settling the 
mentioned sign issues and motivating further theoretical 
investigation. 

We consider a macrospin model corresponding to a 
rectangular MTJ [Fig. [lja)] with nominal free layer di- 
mensions 100x300x2.5 nm 3 , saturation magnetization 
47rMs=4.4 kG, easy axis anisotropy field 2K/Mg = 



2 



Hk=80 Oe, dipole field coupling to the reference layer 
Hjj= 28 Oe favoring the P state, and with parallel 
(P) and antiparallel (AP) resistances i?p=286 f2 and 
Rap =364 £1 respectively, as quoted in Ref. [lij]- The 
magnetization orientation of the free layer, given by the 
unit vector M, is described by the angles and <p de- 
fined in Fig. [Hb) . In the absence of external fields and 
voltage and ignoring any magnetic coupling between the 
free and pinned layers, the equilibrium positions of M 
lie on the easy axis. A schematic of the energy density 
profile at tp = 7r/2, i.e. for in-plane M, is represented in 
Fig. [ljc). We generate a Maxwell-Boltzmann distribu- 
tion of 1000 initial orientations of M corresponding to a 
temperature of 300 K. The time evolution of the ensemble 
of trajectories is found by solving the stochastic Landau- 
Lifshitz- Gilbert (LLG) equation in the presence of both 
Slonczewski's [l[ and field-like [|| spin torque terms, 



T = ajM x (M x nip ) + bjM x m P . 



(1) 



for a macrospin nanomagnet with energy density 

u = -K{M-e\\f-M s {H D +H\\)M.-e\\+K P M-e n . (2) 

The effects of thermal fluctuations during the evolution of 
M are included by using thermal random fields [Tt], EH • 
We describe the angle dependence of the MTJ resistance 
by R(6) = R P + (R A p - Rp) sin 2 0/2. Although the 
results shown here are for an angle independent effi- 
ciency, g(9,p) = r), similar results were obtained by in- 
cluding Slonczewski's MTJ efficiency, g{0,p) = p/(2 + 
2p 2 cos(0)) where p is the spin polarization of the 
tunneling electrons. As proposed by Li et al. Q, we use 
bj = e\V\aj, where V is the voltage and e controls the 
relative amplitude of the two spin torque terms. Ref. [|| 
estimates e ~1 V -1 for typical magnetic materials. We 
model the voltage waveform as a step with 55 ps rise- 
time, v By appropriately choosing the values of e, rj, 
and the damping, a, it is possible to obtain a situation 
where (i) the field-like spin torque is small enough so 
that the P state is still stable and a large majority of 
the thermally distributed initial orientations of M are 
within the stability region of P, as shown in Fig. [ljd), 
and (ii) the Slonczewski's spin torque is too weak to in- 
duce magnetization reversal by itself. Under these condi- 
tions, it is the combination of Slonczewski's spin-torque 
together with thermal fluctuations which eventually push 
the magnetization over the barrier, which has been low- 
ered by the field-like spin torque. From the values of 
Rp and Rap we obtain a zero bias tunneling magnetore- 
sistance TMR«27%, from which the spin polarization 
p «34% and efficiency rj «0.15 can be obtained [2(j. As 
the voltage across the MTJ is increased, the efficiency 
will decrease, possibly by more than 70% when V=l V 
(the voltage used in Ref. [Hj]). 

We first find the time evolution of M for r;=0.15, 
a=0.02, e=1.7 V" 1 , and V=l V. We generate 1000 tra- 
jectories such as the one shown in Fig. [U^a). Then we 
randomly pick ten of them, one by one, with the restric- 
tion that any picked trajectory has to have its switching 




FIG. 2: (Color online) (a) Normalized change in resistance as 
a function of time for a randomly chosen trajectory, (b) Same 
for a set of trajectories with switching times ts distributed 
between and 10 ns. (c) Average of all resistance traces with 
ts > ts after aligning their switching time with the average 
value ts~2.75 ns. Inset: switching via Slonczewski's spin 
torque only (e — 0). All other parameters are kept constant. 



time clearly spaced from all previously picked trajecto- 
ries. This process ensures that the set is representative 
of the whole ensemble. The time dependence of the re- 
sistance for the selected trajectories is shown in Fig. [2jb) 
(compare with Fig. 4 in Ref. We observe (i) a 

slow initial increase in the resistance, (ii) random pre- 
switching fluctuations and reproducible post-switching 
ringing, and (hi) similarity in the behavior of M tra- 
jectories shifted so as to align their switching times. Fur- 
thermore, even though in some traces the transition from 
a slow resistance increase to a fast switch is subtle, clear 
transitions at about 0.2AR max are observable in many 
of the traces. We note that this value corresponds to 
the resistance at the top of the energy barrier (Fig. Did)) 
separating the P and AP states. The average switch- 
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FIG. 3: (a), (b) Average of resistance traces after aligning 
their switching times with ts- In (a) £s=2.7 ns, while in (b) 
?i=7.1 ns. 



ing time for the ensemble of trajectories is t$=2.8 ns, 
close to the value measured in Ref. 

Q at 1.1 V (t s =2.5 
ns), but smaller than what they observed at 1 V (ts=5.3 
ns). However, since the average switching time depends 
in part on the function g(6,p) and the value of e used, 
which are not well known, complete numerical agreement 
with the experimental results is not expected. To confirm 
the randomness of the pre-switching fluctuations and co- 
herence of post-switching oscillations we performed the 
following averaging. From the full sample of 1000 trajec- 
tories we picked the ones with switching time ts in the 
interval ts ± Sts with 5ts being the standard deviation. 
We then shifted all resistance traces to the same switch- 
ing time ts and averaged them. The result (Fig. [2fc)) 
shows that the post-switching ringing is indeed preserved 
but the random fluctuations that precede magnetization 
reversal are averaged out. It also shows that the tran- 
sition from a slow to a sharp resistance increase, which 
occurs close to 0.2AR, is preserved. The inset to Fig.[2Jc) 
allows comparison of this result with the typical switch- 
ing process induced by Slonczewski's spin torque alone 
(e=0). In the latter case one observes a clear build up 
of the precession amplitude proceeding magnetization re- 
versal, while the post-switching oscillations are reduced. 

We tested the sensitivity of our observations to changes 
in the parameters 77, a, and e. Figure [3] shows averaged 
resistance traces for two sets of values of the parameters 
•q, a, and e, illustrating that a wide range of these pa- 
rameters lead to the same general behavior. As long as 
the efficiency is small enough so that Slonczewski's spin 
torque is comparable to the random thermal torques (this 
condition relates r\ and a) , the behavior of the resistance 
is well described by Fig. Even in the case where 

P becomes an energy maximum (for large values of e 
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FIG. 4: Average switching time as a function of applied 
field for voltages between 0.8 and 1.2 V in 0.1 V increments. 
The horizontal dashed line shows the negative field saturation 
value, while the tilted dashed lines are guides to the eye. 



and 77) we observe similar behavior, although the average 
switching times are well below 1 ns, and the probability 
of having an initial orientation of M with a switching 
time of more than 1 ns is negligibly small. However, as 
long as e is chosen so that the P-AP energy barrier is a 
few times larger than the thermal energy, it is possible to 
obtain a wide range of switching times, even exceeding 
the value of 5.3 ns measured in Ref. [14] [Fig.^b)]. We 
observe that as the average switching times increase due 
to a larger P-AP barrier, the amplitude of the random 
pre-switching fluctuations increases, since larger fluctu- 
ations in 9 (and therefore in resistance) are required to 
overcome the barrier. As shown in Fig. ^b) some pre- 
switching oscillations still remain after averaging, but 
their amplitude is much smaller than the amplitude of 
the post-switching oscillations. Therefore, there exists a 
large range of reasonable values for the parameters 77, a, 
and e, which satisfy the conditions described above and 
reproduce the main experimental observations. 

In the present device geometry an application of exter- 
nal magnetic field H\ \ along the magnetic easy axis effec- 
tively modifies the field-like torque constant bj — > bj—Hu 
leaving the aj unchanged. Experiments with separate 
control of the MTJ voltage V and applied field Hu can 
provide valuable information on the spin-transfer torque 
parameters. The behavior of the average switching time 
predicted from simulations is shown in Fig. 2J We used 
the same parameters as in Fig. [5J but varied the voltage 
between 0.8 and 1.2 V in steps of 0.1 V, and applied mag- 
netic fields between -44 and 28 Oe along the easy axis. 
Positive fields favor the P state. The apparent satura- 
tion of the switching time for the 0.8 V and 0.9 V traces 
occurs since we mimic the averaging performed in an ex- 
periment with a 10 ns current pulse (switching events 
with ts > 10 ns are excluded). To understand these sim- 
ulations we note that both applied field and voltage affect 
the energy barrier, with bj decreasing and Hii increasing 



4 



it. Since aj — (hr]/2evM s )V/R where v is the free layer 
volume, and R is the effective device resistance, we get 
bj = (3V 2 with (3 = 2eRvMs/hi]£. The barrier is com- 
pletely eliminated at Hk + Ho +H\ i — b j = which gives 
a crossover field 

if 1 1* = fiV 2 — Hk — Ho, 

(marked by arrows in Fig. |4|) separating the ultrafast 
and normal switching regimes. Although the crossover 
is broad due to finite temperature, one would be able 
to extract H^ from experimental data by extrapolating 
the dependencies ts{H^) above and below the crossover. 
Then plotting Hiu as a function of V 2 one can find 
the coefficient (3 and extract the value of e, i.e., the 
relative strength of the field-like spin torque term. To 
demonstrate the feasibility of such data analysis we fit- 
ted the simulated data and obtained e rj 1.75 V -1 , 
Hk+Hd k, 107.5 Oe for a dataset generated with e = 1.7 
V" 1 and H K + H D = 108 Oe. This method should allow 
verification of the V 2 dependence of the field-like torque 
term and the measurement of its magnitude, character- 



ized in our model by e. 

In conclusion, we show that by considering an addi- 
tional field-like spin torque of magnitude similar to Slon- 
czewski's spin torque, a macrospin description of MTJ 
switching reproduces the "incubation delay" , the lack of 
coherent magnetization precession before reversal, and 
the post-switching oscillations observed in Ref. [lij ■ Our 
analysis suggests an additional experimental procedure 
to measure the voltage dependence of the field-like spin 
torque term: an external magnetic field can be applied 
along the easy axis to cancel or reinforce the field-like 
spin torque, while the average switching time in response 
to voltage pulses of different amplitudes is measured. In 
contrast to previous studies based on careful analysis of 
the antisymmetric component of ferromagnetic resonance 
spectra, or fits of the critical switching voltage to the spin 
torque model, our proposed measurements can provide a 
more direct access to the field-like spin torque term, clar- 
ify its origin, and obtain its voltage dependence in a wider 
range. 
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